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Abstract — Rabi oscillations have been observed in many super- 
conducting devices, and represent prototypical logic operations 
for quantum bits (qubits) in a quantum computer. We use a 
three-level multiphoton analysis to understand the behavior of the 
superconducting phase qubit (current-biased Josephson junction) 
at high microwave drive power. Analytical and numerical results 
for the ac Stark shift, single-photon Rabi frequency, and two- 
photon Rabi frequency are compared to measurements made 
on a dc SQUID phase qubit with Nb/A10a;/Nb tunnel junctions. 
Good agreement is found between theory and experiment. 
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I. Introduction 

SUPERCONDUCTING circuits utilizing Josephson junc- 
tions are a promising technology for quantum compu- 
tation. For the superconducting phase qubit, based on a 
single current- or flux-biased Josephson junction [1], the 
spectroscopy of energy levels was first found some twenty 
years ago [2], while Rabi oscillations [3], [4], multiphoton 
transitions [5], capacitive coupling [6], and coupled oscilla- 
tions [7] have only been demonstrated in the last five years. 
While decoherence remains the dominant obstacle for these 
devices, a number of fundamental physical effects remain to 
be studied. 

In this paper, we explore the physics of the phase qubit 
when irradiated by strong microwave fields. A number of 
related experiments and calculations have appeared in the 
literature. Multiphoton Rabi oscillations of a superconducting 
qubit were first performed in a charge qubit device [8], 
while the ac Stark shift was recently measured in a qubit- 
cavity system [9]. High-power multiphoton spectroscopy of 
a three-junction flux qubit has been accomplished (with up 
to 20 photons) [10]. The importance of diagonal terms in 
the microwave drive was theoretically explored for a SQUID 
qubit in [11], while many multiphoton and multilevel effects 
were recently studied with respect to EIT (electromagnetically- 
induced transparency) [12]. Three-level effects in the phase 
qubit have also been studied numerically [13], perturbatively 
[14], and analytically [15]. 
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Here we compare analytical and numerical results from a 
three-level multiphoton analysis to Rabi oscillation experi- 
ments performed on a dc SQUID phase qubit. We develop 
analytical results to describe the basic Rabi frequencies (the 
oscillation frequency between energy levels) and ac Stark 
shifts (the shift of the optimal resonance condition with 
applied power) of the device, which are in good agreement 
with experiment. In addition to confirming the physics of the 
superconducting phase qubit, these results are important for 
situations where fast qubit operations are performed and for 
high-precision applications. 

II. Model 

The Josephson phase qubit is simply a hysteretic current- 
biased Josephson junction that is well-isolated from dissipation 
in its leads [1]. The Hamiltonian is 



H 



£j(cos7+(//7 c ) 7 ) 



(1) 



where 7 is the gauge-invariant phase difference, p 7 its con- 
jugate momentum, I c the junction's critical current, C its 
capacitance, Ej = hl c /2e, and E c = e 2 /2C. Here we con- 
sider a time-dependent bias current / = J<&, — lac cos u>t. Key 
parameters for the system are the classical plasma frequency 
huj = y/8E c Ej (l - (I dc / 1 c) 2 ) 1 ' ' 4 and the barrier height (of 
the washboard potential) in units of Tiujq, denoted by N s : 



N s = 
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where the approximation is valid for I c — Id c <C I c 
three-level basis, the Hamiltonian has the matrix form 



(2) 



In a 
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A cos tot 




Xoi X 02 
X\\ X12 
X\2 X 2 2 



(3) 

where we have ignored any overall energy shifts. Using 
perturbation theory [16], the energies are approximately 



En = fiUJr, 



1 11 2 465 4 
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with A = (54/Vs) 1 I 2 . The (dimensionless) position matrix 
elements are, to order A 3 : 

a;oo= 3A/2, x i= 2- 1 / 2 (l + llA 2 /4), 

x u = 9A/2, x 12 = 1 + llA 2 /2, (7) 

x 02 = -2- 1 / 2 A, x 22 = 15A/2. 

Finally, the drive amplitude has the explicit form 

1/2 



A= I*c Ej f8Ec 



du) lac 

-^3A' (8) 



and is linear in the microwave current I ac . 

III. Multiphoton Rotating Wave Approximation 

To describe multiphoton processes (see, e.g. [10], [11]) we 
transform the Schrodinger equation to an interaction picture 
using the time-dependent diagonal operator 

ho \ / x o 

ff (t) = | hi +Acoswt xu 
h 2 ) \ x 22 

(9) 

where h — E , hi = E + nhui, h 2 = E Q + 2nfkv, and n 
is the number of photons involved in the — > 1 transition. In 
the interaction-picture the new Hamiltonian is 

H(t) = e ir o H ^ ds l h H{t)e- i ^ Ho{s)ds,h - H (t). 

By using the Bessel-function identity 
12m=-oc eml8 Jrn(z) and averaging H(t) over one period 
T = 2-k/uj of the drive one obtains the multiphoton rotating 
wave Hamiltonian H^: 
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IV. Exact three-level analysis 

From the effective three-level Hamiltonian we can find exact 
solutions to the eigenvalue problem H^ip = Eip. Note that 
in this section we drop the (n) dependence on the matrix 
elements. Letting E = hx, the exact eigevalues of H (n *> are 
given by the roots of the cubic equation: 

x 3 - ( £l + e 2 )x 2 + (eie 2 - + ^?2 + ^02)) x s 

+ l( e 1^02 + £ 2^01 — ^01^12^02) = 0. 

This has the exact solution (see also [15]): 

x n = \(ei + e 2 ) + 2p^ 2 cos(d + 2mr/3) (16) 



(n = 0, 1,2) with 



* = r rccos (~^b 



(17) 



P = + e 2 ) 2 - Jei62 + ^(^oi + ^12 + "02). (18) 
and 

q = -4(ei + e 2 ) 3 + 3 (ei^o 2 + £ 2^oi - ^01^02^12) 
+|(e 1 +e 2 )[e 1 e2-i(0 2 1 +0 2 2 + n 2 2 )]. 

(19) 

The eigenvectors are 

(2(ei - a;„)Oo2 - ^01^12 
-2a;„Oi2 - iloi^o2 
+ 4(ei - x„)x„ 

with normalization 

n = [{2{ei-x n )n m -n m n 12 f 

+ (-2.X„Oi2 - £l i^o 2 ) 2 

-1/2 



(20) 



(21) 
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Using these eigenvectors one can calculate the time-dependent 
amplitudes for a given initial state. That is, if 

- a o (t)\0) + oi(t)|l) + 02(t)|2) (22) 

then, given a„(0), the final amplitudes are 

2 

a n (t)= e-^ t (n|« £ )(^|m)a m (0). (23) 

£,m=0 

The eigenvalues xq, Xi, and x 2 can be used to identify the 
various transitions, and the effective Rabi frequencies are given 
by the differences, e.g. ^r.oi = x — x 2 , fi_R,02 = x 2 — x\. 

V. Phase Qubit Approximations 

Here we describe approximations to the exact eigenvalues 
for the phase qubit, when n = 1. First, for typical drive 
amplitudes the Bessel functions can be approximated by 
Jn{z) ~ (z/2) n /n\. Then, using the position matrix elements 
of the cubic one finds 



4V w h^Axoi, 



n 



(i) 
12 



and 



n 



W ~ Sv^A 2 ^!. 



02 



(24) 



(25) 



We see that is much smaller than Q$ by two small 
factors: A 2 and /u>. Thus, this matrix element can typically 
be neglected. 

In this three-level approach, about each transition there 
will be an avoided crossing of two of the eigenvalues. These 
avoided crossings indicate that oscillations between the origi- 
nal energy levels will occur, and with a frequency given by the 
eigenvalue difference. This frequency can be approximated by 
using the characteristic Rabi formula such as 



Qr.,01 = [^a,oi + ( w oi + Awqi - uj) z 



(26) 



where Aw 01 is a possible shift of the resonance position due 
to the third level. 
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ACJQI 



(27) 



The effect of this non-resonant level can be found perturba- 
tively by either adiabatic elimination [12] or by a systematic 
expansion of the exact eigenvalues [15]. Here we simply state 
the final results. For the — > 1 transition the minimum Rabi 
frequency occurs for 

2(wqi - W12) ' 

thus identifying Acjoi = lo — uJqx- This shift of the optimal 
transition point is known as the ac Stark shift in atomic 
physics, and has the characteristic dependence f2 2 /(Ao;), 
where Aui is the frequency difference of the off-resonant 
transition (here lj%2) and the drive (here w woi) and 57 the 
characteristic matrix element. In addition, the minimal Rabi 
frequency is actually decreased from its nominal value of fl^ : 



s 'fi,01 



n 



(i) 
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[n, 



(1)12 



4(w i 



W12) 



(28) 



The two-photon 



2 transition, when can be 



neglected, is produced by two off -resonant single-photon tran 
sitions — > 1 — > 2, with a Rabi frequency of 



(i) 

R,02 



(cjoi - W12) 



(29) 



In addition, a similar Stark shift is found in the 
transition, with the resonance condition: 



W02/2 



(30) 



4(u> i -W12) 

Note that in d27l i and d30b the transition frequencies appear 
on both sides of the equation, and thus should be determined 
self-consistently; this is typically a small effect. 

For these expressions to be valid, the condition O5Y <C 
ujqi— uo 12 must be met. While this is not the case for the highest 
power results of the following section, these perturbative 
results do indicate the main physics. In the following, we will 
compare both these and the full three-level multiphoton effects 
predicted by (TToT l with experiment. 

VI. Experiment 

Measurements were made on a dc SQUID phase qubit 
with two Nb/AlOaYNb Josephson junctions in a 6-turn 6 
nH loop. The devices were fabricated at Hypres, Inc., using 
a trilayer process Q and the qubit junction had an area of 
(10/im) 2 . The SQUID was dynamically initialized to one of its 
metastable wells and simultaneously current- and flux-biased 
in such a way to study the properties of one of its junctions; 
further experimental details can be found in [17]. Microwave 
spectroscopy [2] performed on the device yielded single and 
two-photon transitions in quantitative agreement with quantum 
mechanical spectra of a single junction with parameters C = 
4.46 pF and I c — 17.821 fiA; a typical scan in frequency 
and bias current is shown in Fig. 1. Applying microwaves of 

1 Mention of commercial products or services in this paper does not imply 
approval or endorsement by NIST, nor does it imply that such products or 
services are necessarily the best available for the purpose. 
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Fig. 1. Experimental microwave spectroscopy of a Josephson phase qubit, 
scanned in frequency (vertical) and bias current (horizontal). Dark points 
indicate experimental microwave enhancement of the tunneling escape rate, 
while white dashed lines are quantum mechanical calculations of (from right 
to left) u> i, ^02/2, u>i2, ^13/2, and lu 2 3 ■ 




Fig. 2. Rabi frequency f2ij,oi of the one-photon -t 1 transition as function 
of microwave current I ac - The dots are experimental data, the solid line 
predictions from the three-level model, and the dashed lines are the lowest- 
order results 1241 (top) and second-order (28) (bottom) perturbative results. 
The inset shows Rabi oscillations of the escape rate for I ac = 16.5 nA. 
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Fig. 3. The ac Stark shift Ao;oi of the one-photon — » 1 transition as 
function of microwave current I ac . The dots are experimental data, the solid 
line predictions from the three-level model, and the dashed line perturbative 
results. The inset shows the oscillation frequency Or qi as a function of the 
level spacing u>oi for I ac = 5.87 nA and the fit using {26} to obtain Qr qi 
and AaJrji. 
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Fig. 4. Rabi frequency of the two-photon — > 2 transition as 

function of microwave current I ac . The dots are experimental data, the solid 
line predictions from the three-level model, and the dashed line perturbative 
results. Inset shows Rabi oscillations of the escape rate for I ac = 16.5 nA. 



constant frequency oj/2it = 5.9 GHz, resonances were seen 
for bias currents matching a single or two-photon transition. 
Near Id c = 17.660 /iA the single-photon — > 1 transition was 
observed, corresponding to energy level spacings of luqi/2tt = 
5.9 GHz and uj^jl-K = 4.99 GHz. For the same microwave 
frequency but at a smaller bias Id c — 17.636 /iA the two- 
photon 0^2 transition was observed, with levels spacings 
of cj i/2tt = 6.23 GHz and uj 12 /2tt = 5.57 GHz. 

About each transition, and for several microwave powers, 
Rabi oscillations were observed in the tunneling rate of the 
junction from the zero to finite voltage state [4] (see inset 
to Fig. 2). The Rabi oscillation frequency, shown in Fig. 2, 
was found by varying the bias current and noting the level 
spacing at which maximal oscillations of the escape rate and 
minimal Rabi frequency was located (see inset to Fig. 3). Also 
shown are the exact three-level results and approximate results 
of the previous section. To produce these curves a single 
fitting parameter was used to calibrate the applied microwave 
current I ac as a function of the applied microwave power 
P. This calibration (specifically, I ac = 58.7nA(P/mW) 1 / 2 ) 
was then used for the horizontal axis. Note also that the 
numerical values for a>oi and u>\2 were used in the theory, as 
the perturbative results ©-(O are accurate only for N s > 3 
[16] while here N s w 2.21. 

Using this calibration, the ac Stark shift Aujqi = u> — woi 
(derived from the bias current) is plotted in Fig. 3. The solid 
curve is the three-level prediction, with no fitting parameters, 
while the dashed curve is the perturbative result d27t . We see 
good agreement with the three-level results, even for powers 
where the perturbative result fails. 

Finally, the Rabi frequency of the two-photon 0^2 transi- 
tion is shown in Fig. 4. Here the solid curve is the three-level 
prediction while the dashed curve is the approximation d29l ). 
For small powers, the experimental results are likely domi- 
nated by systematic errors in the tunneling rate measurements. 
For large powers we find that while the perturbative result 
fails, the full three-level prediction is in excellent quantitative 
agreement with the data. 



VII. Conclusion 

In conclusion, we have examined the behavior of the su- 
perconducting phase qubit when irradiated with microwaves 
at high power. We have identified two effects that are new 
to this qubit: multiphoton Rabi oscillations and the ac Stark 
shifts. We have also shown that a simple three-level theory can 
explain, with minimal fitting parameters, both of these subtle 
experimental effects. This is quite remarkable, as dissipation, 
decoherence, and tunneling are all signficant in the same 
parameter regime for these devices. Finally, these results 
indicate that high-power Rabi oscillations can in principle be 
used to characterize all of the relevant matrix elements needed 
for high-fidelity microwave control of superconducting phase 
qubits. 
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